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Abstract 

We compute the form factors for the B —)• Kl^l~ semileptonic decay process in lattice QCD 
using gauge-field ensembles with 2+1 flavors of sea quark, generated by the MILC Collaboration. 
The ensembles span lattice spacings from 0.12 to 0.045 fm and have multiple sea-quark masses 
to help control the chiral extrapolation. The asqtad improved staggered action is used for the 
light valence and sea quarks, and the clover action with the Fermilab interpretation is used for 
the heavy b quark. We present results for the form factors /+(g^), /o((?^), and where is 

the momentum transfer, together with a comprehensive examination of systematic errors. Lattice 
QCD determines the form factors for a limited range of and we use the model-independent 
z expansion to cover the whole kinematically allowed range. We present our final form-factor 
results as coefficients of the z expansion and the correlations between them, where the errors on 
the coefficients include statistical and all systematic uncertainties. We use this complete description 
of the form factors to test QCD predictions of the form factors at high and low q^. 
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i Present address: Laboratory of Biological Modeling, NIDDK, NIH, Bethesda, Maryland, USA 
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I. INTRODUCTION 


Flavor-changing nentral-cnrrent interactions (FCNC) place important constraints on 
physics beyond the Standard Model. In the Standard Model, tree-level FCNC contribntions 
vanish by the Glashow-Iliopolons-Maiani (GIM) mechanism. Even at the one-loop level, 
the GIM mechanism snppresses these amplitudes, as do factors of the Gabibbo-Kobayashi- 
Maskawa (GKM) mixing matrix. Thus, new-physics effects may be substantially larger than 
the small Standard-Model contribution and, hence, observable. In this paper, we present an 
unquenched lattice-QGD calculation of the amplitudes for the FGNG process B —)■ Kl^l~. 
Within the Standard Model and beyond, three form factors can arise, and we present results 
for all three. This work is part of a larger program by the Fermilab Lattice and MILG 
collaborations to calculate form factors for exclusive semileptonic B decays needed to test 
the Standard Model and search for new physics, all of which use the same lattice actions 
and parameters. It builds upon our previous work on charged-current semileptonic B de¬ 
cays, B —)■ ttIv 0,0 and B —)■ 0-0, which are used to determine the GKM matrix 

elements |I4fe| and \Vcb\ 0 - If is also part of a suite of form factors needed for searching for 
new physics in rare semileptonic 5-decay processes such as 5 —)■ 0 , 5 —)- Dtu 0 

and Bg —?• 0. 


Experimental research on rare 5-meson decays is active |10l. Illl|. The BaBar, Belle, and 


GDF collaborations have measured the differential branching ratio, the forward-backward 
asymmetry and other observables for both 5 —)■ Kl^l~ and 5 —)■ K*l^l~ decays |I2^ 


16 


The LHGb Gollaboration has reported more precise results for the B^ 

decays [ 13-0 


5 ± ^ K^+l 
the future 20 


—)■ K^l^l and 
The high-intensity 5 factories will also have results in 


and 


Thus, it is timely to improve the precision of the theoretical calculation 
of these processes. Recently, the HPQGD Gollaboration published the hrst three-flavor 
lattice calculation of 5 —)■ Kl^l~ j2l|, also analyzing the phenomenological implications 22 
Three-flavor results for the modes with vector mesons in the hnal state, 5 —)■ 

Bg —>■ have also been presented [0 |. 

The theoretical description of the 5 —)■ process is based on the operator-product 

expansion, which leads to a low-energy effective Hamiltonian [24-23]. Amplitudes are ex¬ 
pressed in terms of Wilson coefficients that encode the high-scale physics and hadronic 
matrix elements that capture the low-energy nonperturbative QGD contributions. Hadronic 
matrix elements of local operators can be parameterized in terms of form factors. The un¬ 
certainties in the form factors are an important source of error in the theoretical predictions 
of the observables mentioned above. In order to calculate the form factors, one requires 
knowledge of nonperturbative QGD dynamics, and lattice QGD is the tool of choice. We 
focus on 5 —)■ Kl^l~, rather than 5 —)• because the former is a “gold-plated” decay 

with a stable hadron (under strong interactions) in the hnal state. In the vector-meson case, 
the K* is unstable, and the resonance would have to be distinguished from Ktt states. 

The goal of this work is to improve our knowledge of the 5 —)■ Kl^l~ form factors. We use 
the three-havor gauge-held ensembles generated by the MILG Gollaboration with dynamical 
up, down, and strange quarks. We extrapolate our lattice simulation data to the physical 
light-quark masses and continuum using SU(2) chiral perturbation theory formulated for the 
process 5 —)■ Kl^l~. Because the strange-quark mass is integrated out of the SU(2) theory, 
the use of SU(2) xPT, rather than SU(3), improves the convergence of the chiral expansion, 
thereby reducing the systematic uncertainty due to the chiral-continuum extrapolation. On 
currently available lattices, we directly obtain the form factors at large momentum transfer 
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(low recoil), > 17 GeV^. Following Ref. [I|, we use the ^ expansion to extend the 
lattice-QCD calculation to the full range of Compared with the work of the HPQCD 
Collaboration 2l|, l22|, we use twice as many ensembles, covering a larger range of lattice 


spacings and using lighter sea-quark masses. In particular, the smallest lattice spacing and 
the smallest light-quark mass in our ensemble set are smaller by a factor of two compared to 
the set used by HPQCD. In addition, we use the Sheikholeslami-Wohlert (SW) action [ 2 ^ 
with the Fermilab interpretation 12911 for the 6-quark, while the HPQCD Collaboration uses 
a nonrelativistic QCD b quark As discussed below, details of the chiral-continuum 

extrapolation and the 2 ; expansion also differ. 

This paper is organized as follows. In Sec. HIl we define the form factors for the B —)■ 
Kl^l~ decay. We then describe the lattice ensembles used in our simulations. We also 
discuss the formalism for the light and heavy quarks. In Sec IHIl we present the numerical 
analysis. This section has four parts. We first present results for the simulated K and 
B meson masses. Next, we determine the lattice form factors from two-point and three- 
point correlation functions. We then correct our form-factor data to account for the slight 
difference between the simulated 6-quark mass and the physical 6-quark mass. Last, we 
extrapolate the lattice simulation results to the chiral and continuum limits with SU(2) 
heavy meson rooted staggered chiral perturbation theory (HMrSyPT). In Sec. IIVI we analyze 
the systematic errors in our calculation and give a complete error budget for the range of 
> 17 GeV^ accessible in our numerical simulations. In Sec. we extrapolate 


We present our 


momenta g 

our form factors from low to high recoil using the z expansion j3ll-34|. 
hnal results for /+(g^), /o(g^), and /T(g^), including statistical and all systematic errors, as 
coefficients of the expansion and the correlations between them; this provides a complete 
description of our form factors valid over the entire kinematic range. In Sec. I VII we use 
these results to test predictions for the form factor from heavy-quark symmetry at high g^ 
and from QCD factorization at low g^. Finally, we compare our form factors with other 
lattice-QCD and light-cone-sum-rule results, and present an outlook for future work, in 
Sec. Em 

Preliminary results have been reported in Refs. j35-37|. Here we present a full analysis 


that includes the tensor-current form factor and complete systematic error budgets. 


II. LATTICE CALCULATION 

In this section, we present the methods and ingredients used in this work. We give the 
definitions of the form factors for the B —)■ Kl^l~ process and their relation to physical 
observables in Sec. Ill Al We then describe the lattice actions and parameters used for gluon 
and fermion helds in our simulation in Sec. Ill HI Finally, we define the lattice currents in 
Sec.Hm 


A. Matrix elements and form factors 

An operator-product expansion analysis of the B —)■ Kl^l~ decay in the Standard Model 
shows that two currents, a vector current = 57^6 and a tensor current = isa^^'b, 
contribute to the 6 —)■ s process at lowest order [l^. In general Standard Model extensions, 
a scalar current S = sb can also arise. The matrix elements of the vector, tensor, and 
scalar current are characterized by three form factors /+(g^), /o(g^), and /r(g^), which are 
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defined via 


/ A/f2 _ /\/f2 \ ]\/f2 _ T\/f2 

{K\srh\B) = U{q^) + F - + /o(g^) ^ ^ q^ (2.1) 

(A1iS(r'“'6|B) = - p-A.-), (2.2) 

(A'|56|B) = (2,3) 

mb - rris 

where p and k are the 5-meson and kaon momenta, respectively, and q = p — k is the 
momentum carried off by the leptons. The Ward identity relating the matrix element of 
a vector current to that of the corresponding scalar current ensures that /o is the same in 
Eqs. fl2.ip and fl2.3p . 

For the analysis that follows, it is convenient to write the vector-current matrix element as 


{K\srb\B) = [v^fiiiEK) + k^J^iEK)] , (2.4) 


where = p^/M b is the four-velocity of the B meson, k^ = k^ — {k ■ v)v^, and Ek = 
V ■ k is the kaon energy in the 5-meson rest frame. From energy-momentum conservation, 
q^ = — 2 MbEk- We obtain /||(5/f) and /±(5i^) from the temporal and spatial 

components of the matrix element of the vector current: 


MEk) 

f±{EK) 


{K\sy%\B) 

{K\sy%\B) 

y/2M^k^ 


(2.5) 

( 2 . 6 ) 


Similarly, we obtain the tensor form factor f-r from 


Mq^) 


Mb + Mk {K\sa^^b\B) 
^/2Mb ^/2MbE 


(2.7) 


Finally, the vector and scalar form factors /+ and /o can be obtained from 


fM) = :7= [/ii(Ak) + {Mb - EK)fL{EK)\ , (2.8) 

Mf) = - EK)h{EK) + (El - MI)!b{Ek)] . (2.9) 


Equations fl2.8p - fl2.9p satisfy the kinematic constraint, /+(0) = /o(0), automatically. At low 
recoil, the form factor f± gives the dominant contribution to /+. 

Physical observables can be described in terms of the form factors, if we neglect non- 
factorizable contributions. For example, the Standard-Model differential decay rate for 


5 ^ Kl+l- is [24, |38|, [39 


dT 

dq^ 


Gla^ 


\V.„V, 


ts I 




\k\/3, 




\CfoU{q^ 






\Cfofo{q^) 



CfUiq^) + 2Cf 


rUb + rus 
Mb + Mx 



( 2 . 10 ) 
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TABLE I. Parameters of the QCD gauge-field ensembles and light valence-quark masses used in 
this work, lattice spacing a, lattice size Ng x Nt, sea-quark masses am'i and light-valence 

mass ami, daughter mass amh, the number of configurations and sources denoted as A"conf x Agro 
and the box size times the pion mass. On all ensembles but one, we use the same light valence- 
and sea-quark mass. (The only exception is on the a ~ 0.09 fm ensemble with = 0.0465, where 
the light valence-quark mass is 0.0047 instead of 0.00465.) On the a ~ 0.12 fm and a ~ 0.09 fm 
ensembles we also use the same valence and sea strange-quark mass. On the a ~ 0.06 fm and 
a ~ 0.045 fm ensembles, we use slightly different valence strange-quark masses than in the sea; the 
valence masses are tuned to be closer to the physical value. The values of are taken from 

Refs. 0, 52|. The gauge-held configurations can be downloaded using the digital object identiher 
(DOI) links provided in Refs. 66-81|. 


Rs a (fm) 

N^xNt 

am'i 

am'h 

ami 

amh 

Aconf X Agj-c 

M^L 

0.12 [66,^ 

20^x64 

0.01 

0.05 

0.01 

0.05 

2259x4 

4.5 

0.12 

20^x64 

0.007 

0.05 

0.007 

0.05 

2110x4 

3.8 

0.12 

24^x64 

0.005 

0.05 

0.005 

0.05 

2099x4 

3.8 

0.09 [70-7^ 

28^x96 

0.0062 

0.031 

0.0062 

0.031 

1931x4 

4.1 

0.09 [7^ 

323x96 

0.00465 

0.031 

0.0047 

0.031 

984x4 

4.1 

0.09 [74,1^ 

403x96 

0.0031 

0.031 

0.0031 

0.031 

1015x4 

4.2 

0.09 [7^ 

643x94 

0.00155 

0.031 

0.00155 

0.031 

791x4 

4.8 


0.06 [77, 7^ 
0.06 [79,^ 

483x144 

643x144 

0.0036 

0.0018 

0.0180 

0.0180 

0.0036 

0.0018 

0.0188 

0.0188 

673x4 

827x4 

4.5 

4.3 

0.045 [81J 

643x192 

0.0028 

0.0140 

0.0028 

0.0130 

801x4 

4.6 


where Gp, a, and Vtq are the Fermi constant, the (QED) hne structure constant, and CKM 
matrix elements, respectively, |fc| = is the kaon momentum in the S-meson rest 

frame, and = 1 — Amf/q'^, with mi being the lepton mass. The Cf® are effective Wilson 


coefficients [2^; we follow the notation of Ref. 27(| in Eq. fl2.10p . When corresponds to 


a charmonium resonance, further contributions must be added to Eq. fl2.10p . Beyond the 
Standard Model, the expression can become more complicated, but /+(g^), and /o(<?^) 

still suffice. 


B. Actions and parameters 


Our calculations employ the Nf = 2 + 1 flavor gauge configurations generated by the 

4l| . which include the effects of dynamical u, d, and s quarks. The 
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one-loop improved Liischer-Weisz action is used for the gluon fields, which leads to lattice 
artifacts of O(a.^a^) jd^ . (The gluon-loop correction is included ji^, but not that of the 
quark loop |44|.) 

For light quarks (u, d and s), these configurations employ the tadpole-improved stag¬ 
gered action (asqtad) j45-M|, leading to discretization errors of 0{asa‘^) and [H^. The 

sea quarks are simulated with the fourth root of the staggered fermion determinant. Several 
theoretical and numerical analyses support the idea that this procedure yields continuum 
QCD as the lattice spacing a —)■ 0 |52l-l65j|. 

Table [I] summarizes the properties of the ensembles used in this work. We use the asqtad 
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TABLE II. Parameters used in the simulation of the heavy b quark We list the clover coefficient 

csw, input 6-quark hopping parameter and rotation coefficient di. 


a (fm) 

am'i 

Csw 

< 

di 

0.12 

0.01 

1.531 

0.0901 

0.093340 

0.12 

0.007 

1.530 

0.0901 

0.093320 

0.12 

0.005 

1.530 

0.0901 

0.093320 


0.09 

0.0062 

1.476 

0.0979 

0.096765 

0.09 

0.00465 

1.477 

0.0977 

0.096708 

0.09 

0.0031 

1.478 

0.0976 

0.096690 

0.09 

0.00155 

1.4784 

0.0976 

0.096700 


0.06 

0.0036 

1.4287 

0.1052 

0.096300 

0.06 

0.0018 

1.4298 

0.1052 

0.096300 

0.045 

0.0028 

1.3943 

0.1143 

0.08864 


ensembles at four lattice spacings: a ~ 0.12 fm, a ~ 0.09 fm, a ~ 0.06 fm, and a ~ 0.045 fm. 
The volumes of the lattices are large enough {M^L > 4) to suppress finite-volume effects. 
The strange sea-quark mass is tuned to be close to its physical value. The light-to-strange 
sea-quark mass ratios range from am'Jam'^ = 0.2 down to 0.05, to facilitate reliable chiral 
extrapolations. On the a ~ 0.12 fm and a ~ 0.09 fm ensembles, we use unitary data, with 
the light and strange valence-quark masses set equal to the corresponding sea-quark masses, 
with one exception. On the a ~ 0.06 fm and a ~ 0.045 fm ensembles, however, we use 
valence strange-quark masses that are closer to the physical value and, thus, differ slightly 
from the strange-quark mass in the sea. 


On each conhguration, we compute the correlation functions starting at four different 
source locations, to increase the available statistics. We hrst translate the gauge held by a 
different random four-vector on each conhguration and then hx the spatial source locations 
at a; = 0 and the temporal source locations at f = 0, W/4, Nt/2, and SNt/A. The correlation 
between the results from diherent source locations is weak. The random translation of the 
gauge held reduces autocorrelations between successive conhgurations. 

For the heavy b quark, we use the Sheikholeslami-Wohlert (SW) action 


28 with the 


Fermilab interpretation 1291 . The lattice action and currents are matched to the continuum 
QCD action via HQET |82j. The heavy-quark action can be systematically improved to 
arbitrarily high orders in l/m^, or, equivalently, a, W including higher-dimensional operators 


in the lattice action 


82, 83 and currents [8: 


In this work, we remove the leading 


discretization errors in the action by tuning the hopping parameter k and clover coefficient 
csw- We hx the bare 6-quark mass by tuning the value of Kb to reproduce the spin-averaged 


meson kinetic mass as in Ref. 


We use the tadpole-improved tree-level value for 


Csw = where uq is obtained from the fourth root of the plaquette. We also remove the 
leading discretization error in the vector and tensor currents; see Sec. Ill Cl The values of the 
parameters for 6 quarks used in our simulations are listed in Table [TTl 

To extrapolate the form factors calculated on the lattice to the continuum limit, we need 
a unihed scale to compare the results from diherent spacings and convert to physical units. 
We do so with the scale ri which is dehned such that r^F(ri) = 1.0 87, Here F{r) is 
the force between static quarks at distance r. We hrst determine the relative scale ri/a on 
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TABLE III. Relative scales ri/a used in this work, for corresponding values of /3 [ 1 , 52|. 
statistical and systematic errors on ri/a are both 0.1-0.3% [^. We also list the Goldstone 
mass (M,r) and root-mean-square (RMS) pion mass here. 


« a (fm) 

am'Jam'y^ 

/3 

ri/a 

M^(MeV) 

MpP^^{MeY) 

0.12 

0.01/0.05 

6.760 

2.739 

389 

532 

0.12 

0.007/0.05 

6.760 

2.739 

327 

488 

0.12 

0.005/0.05 

6.760 

2.739 

277 

456 

0.09 

0.0062/0.031 

7.090 

3.789 

354 

413 

0.09 

0.00465/0.031 

7.085 

3.772 

307 

374 

0.09 

0.0031/0.031 

7.080 

3.755 

249 

329 

0.09 

0.00155/0.031 

7.075 

3.738 

177 

277 

0.06 

0.0036/0.018 

7.470 

5.353 

316 

340 

0.06 

0.0018/0.018 

7.460 

5.307 

224 

255 

0.045 

0.0028/0.014 

7.810 

7.208 

324 

331 


The 

pion 


each ensemble, and then interpolate ri/a with a smooth fnnction of the gange conpling /3; 
the smoothed ri valnes are independent of the light sea-qnark mass. (The explicit form of 
the smoothing fnnction is given in Ref. Q.) In this paper, we choose a mass-independent 
scheme for ri/a, so that it is the same for all sea masses with the same approximate lattice 
spacing. We nse the valnes of ri/a to convert all lattice qnantities to ri nnits. We can then 
combine resnlts from different ensembles and perform a chiral-continnnm extrapolation. The 


physical valne ri = 0.3117(22) fm [52, l89| is determined by reqniring that the continnnm 


limit of the pion decay constant at the physical qnark masses takes the PDG valne 
The RBC-UKQCD Collaboration also reported the physical valne ri = 0.323(8)(4) fm in 


Ref. |90|. This resnlt is consistent with the one we use, but less precise. The values of ri/a 


used in this work are provided in Table IIIIl 


C. Definition of currents 


We dehne the current operators on the lattice as in 


Refs. [11,0: 


Vl'ix) = (2.11) 

T^''{x) = T„(a:)cr^^C^^(a;)x(a;), (2.12) 


where the matrix C = x{x) is the one-component staggered fermion 

held. The clover 6-quark held is rotated to remove discretization errors of order a from the 
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>1' = (1 + adi7 ■ Dia,))/!, 


(2.13) 


where ip is the held in the Fermilab action (for the b quark), ZJiat is the symmetric, nearest- 
neighbor, covariant diherence operator, and di is adjusted to remove discretization errors. 
In practice, we set the rotation coefficient di to its tadpole-improved tree-level value: 



Mo 


1 



2 -|- rriQa 2(1 -|- mga) 


(2.14) 



















where moo is the bare lattice &-quark mass. The index ^ in Eqs. fl2.1ip and fl2.12p corresponds 
to taste, and it is contracted with another taste index in the heavy-light operators coupling 
the B meson to the vacuum [^ . 

To calculate the form factors on the lattice, we have to define currents with the correct 


continuum limit. As in earlier work [H, 1^ , we define 

= Zv.V>^, 


(2.15) 

(2.16) 


where {V, T} and {E, T} are the continuum and lattice current operators, respectively. We 
use a mostly nonperturbative renormalization procedure to obtain the Z factors 


Zj - pj^lZy4Zvi, 


(2.17) 


where Zyi^ and Zy4^ are computed nonperturbatively, and the remaining factor pj is calcu¬ 
lated at one-loop order in mean-held improved lattice perturbation theory j^. 

The light-light renormalization factor Zy4^ is calculated nonperturbatively from the charge 
normalization condition of a cs meson: 


Z-i = / d^x{D,\V,Ux)m 


(2.18) 


as in Ref. [891], but with random color wall sources and higher statistics, leading to the 


values listed in Table IIVI The result for Zy4 is insensitive to the mass of the spectator 
quark in the correlation function, so we use a heavy charm quark to improve the statistical 
errors. The heavy-heavy renormalization factor Zy4 is computed analogously from the 
charge normalization condition of the B meson using data generated for our B ^ Dlv 
analysis [^. We compute Zy4^ on the same jackknife samples as the form factors and 
propagate the statistical error directly throughout the remainder of the analysis. The values 
of Zy4 are shown in Table IIVI 

bb I 

The remaining factor pj (here, J = is close to unity |8i 


94 


because most 

of the radiative corrections, particularly those from tadpole diagrams, cancel among the Z 
factors in Eq. fl2.17p . We expand the factor pj perturbatively as 


Pj = 1 + a,(q’)p'f + 0(0^), 


(2.19) 


where is the QCD coupling [95|, |96|. Details of the one-loop perturbative calculation 
will be given in a separate publication j^; the values used here are listed in Table |Vl In 
practice, we evaluate the coupling in the V-scheme 95j, l98j at the scale q* = 2/a in mean- 
field improved lattice perturbation theory. For py4 we find that the one-loop corrections 
are less than 1%, while for pyi they range from 1.5% to 2.6%. The tensor current is scale 
dependent, and we renormalize it at the scale p = rub (where according to the Fermilab 
prescription m 2 = mb). We hnd that for p^ the corrections range from 3% to 6%. 

Because pj is computed separately from the correlation functions, we used it to introduce 
a blinding procedure (as in many B physics experiments) to reduce subjective bias. Those 
of us carrying out the perturbative calculation [^ multiplied pj by a constant prefactor. 
Only after we hnalized the choices made in our analysis, including tests and estimates of 
systematic uncertainties, was the prefactor revealed to the rest of the collaboration and 
removed from the results reported here. 
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TABLE IV. The flavor-conserving renormalization factors Zyt and Zy4 used in this work. Errors 

S'S ^ bb 

shown are statistical. 


Ri a (fm) 

anil 

am'h 


Zy4 

Zy4 

^bb 

0.12 

0.01 

0.05 

0.0901 

1.741(3) 

0.5065(57) 

0.12 

0.007 

0.05 

0.0901 

1.741(3) 

0.5119(75) 

0.12 

0.005 

0.05 

0.0901 

1.741(3) 

0.5026(71) 

0.09 

0.0062 

0.031 

0.0979 

1.777(5) 

0.4482(57) 

0.09 

0.00465 

0.031 

0.0977 

1.776(5) 

0.4694(100) 

0.09 

0.0031 

0.031 

0.0976 

1.776(5) 

0.4608(94) 

0.09 

0.00155 

0.031 

0.0976 

1.776(5) 

0.4491(116) 

0.06 

0.0036 

0.018 

0.1052 

1.808(6) 

0.4196(101) 

0.06 

0.0018 

0.018 

0.1052 

1.807(7) 

0.4100(103) 

0.045 

0.0028 

0.014 

0.1143 

1.841(6) 

0.3564(65) 


TABLE V. Matching factors py4, pyi, and pT calculated at one loop in tadpole-improved lattice 
perturbation theory. Here, pT brings Jt to the MS scheme at p = m 2 , and m 2 should be interpreted 
as the pole mass. 

ss a (fm) 

am'i 

am'h 

amh 

< 

Py4 

dyi 

Pt{p = m 2 ) 

0.12 

0.010 

0.050 

0.050 

0.0901 

1.0071 

0.9737 

1.0334 

0.12 

0.007 

0.050 

0.050 

0.0901 

1.0071 

0.9737 

1.0333 

0.12 

0.005 

0.050 

0.050 

0.0901 

1.0072 

0.9738 

1.0333 

0.09 

0.0062 

0.031 

0.031 

0.0979 

0.9997 

0.9759 

1.0366 

0.09 

0.00465 

0.031 

0.031 

0.0977 

0.9998 

0.9759 

1.0364 

0.09 

0.0031 

0.031 

0.031 

0.0976 

0.9999 

0.9758 

1.0364 

0.09 

0.00155 

0.031 

0.031 

0.0976 

0.9999 

0.9757 

1.0364 

0.06 

0.0036 

0.018 

0.0188 

0.1052 

0.9956 

0.9792 

1.0432 

0.06 

0.0018 

0.018 

0.0188 

0.1052 

0.9956 

0.9792 

1.0433 

0.045 

0.0028 

0.014 

0.013 

0.1143 

0.9943 

0.9843 

1.0588 


III. ANALYSIS 


In this section, we present our form-factor analysis and results. In Sec. IIII Al we obtain 
the i?-meson and kaon masses and energies by htting two-point correlation functions. In 
Sec. IIIIBl we extract the lattice form factors from ratios of three-point over two-point 
correlation functions. In Sec. IIII Cl we slightly shift the full set of lattice form-factor data 
from the simulated to the physical value. In Sec. IIII Dl we carry out the chiral-continuum 
extrapolation by htting the form factors to the expression derived in heavy meson rooted 
staggered chiral perturbation theory (HMrSxPT). 
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A. B and K meson masses 


We extract meson masses and energies from two-point correlation functions defined at 
Euclidean time t: 

k) = Y,{Op{x, f)(9j,(0, (3.1) 

X 

where the subscript P denotes the K or B pseudoscalar meson in the interpolating operator. 
For the kaon we use a local interpolating operator. For the B meson we use the wave fun ction 
for bottomonium given by the Richardson potential model 99] as explained in Refs. [I^ 
13 . We generate correlators with kaon three-momenta k = 27r(0, 0, 0)/L, 27r(l, 0, 0)/L, 


27r(l, 1, 0)/L, and 27r(l, 1, 1)/L. 

The meson masses and energies are extracted from the large-f behavior of the two-point 
correlation functions. By inserting a complete set of states, two-point correlation functions 
can be decomposed into a sum of energy levels as 


4 ' 2B'”' 




(3.2) 


The amplitudes of terms with odd m oscillate in time as (—1)”^(*+^) and are due to opposite- 
parity-state contributions to staggered correlators. Figure [1] shows sample kaon and R-meson 
scaled correlators [C 2 {t) — C^\t)]/C^\t) on the a ~ 0.12 fm ensemble with m'l = O.lm)^ 
and momentum A; = 0, where 



|(0|Op|P(0))p .(0), 


2E 


( 0 ) 


(3.3) 


is the ground-state contribution determined by our £t. The opposite-parity-state contribu¬ 
tion is insignificant for the zero-momentum kaon but is visible for the B meson. We employ 
a simple strategy to fit the two-point correlators because the statistical errors in the kaon 
and R-meson energies contribute little to the errors in form factors, which stem primarily 
from the three-point correlators. For the kaon correlators, we perform two-state hts that 
include the ground state and a same-parity excited state. For the R-meson correlators, we 
perform three-state hts including the ground state, its excited state, and the lowest-lying 
opposite-parity state. 

We use a single-elimination jackknife analysis to estimate the statistical errors in this 
work. We first average the correlation functions generated from the four sources at 0, N^/A, 
Nt/2, and 3Nt/A. We £t C' 2 (t) in an interval t G [tminPmax], taking correlation from time 
slice to time slice into account. In general, we choose fmax so that the fractional error in 
the correlator remains below 4%. We choose fmin such that we obtain a good correlated 
p value. We use the same interval [fmin, ^max] for all kaon or R-meson hts at a given lattice 
spacing, and use similar physical distances for [fmin, ^max] on the four lattice spacings. These 
ht ranges are given in Table IVTl We use a 2-|-l-state ht for the R meson in this paper and 
hnd consistent results with the H-1-state, larger fmin ht of Ref. j^. 

Figure [ 2 ] shows sample R-meson and kaon correlator hts versus fmin for hxed fmax on the 
same a ~ 0.12 fm ensemble as in Fig. [TJ The ht results and errors are stable versus 
and show no evidence of residual excited-state contamination. 

For kaons with nonzero momentum, we can either extract the energy from two-point 
correlation functions with nonzero momentum, or we can use the kaon mass from the zero- 
momentum correlator and the continuum dispersion relation, -|- k^. Figure [3] 
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FIG. 1. Scaled correlator [C 2 {t) — c‘^\t)]/C^\t) as a function of time t on the am'Jam'i^ = 
0.005/0.05 a ~ 0.12 fm ensemble at the unitary point. The oscillating opposite-parity-state con¬ 
tribution is clearly visible in the i?-meson correlator (left), but it is small in the zero-momentum 
kaon correlator (right). 


TABLE VI. Fit ranges [tmin; ^max] used in the kaon and i?-meson mass and energy fits. 


PS a (fm) 

kaon 

B meson 

0.12 

[7,30] 

[3,15] 

0.09 

[10,35] 

[5,20] 

0.06 

[17,60] 

[7,30] 

0.045 

[20,90] 

[8,40] 


shows a comparison of the kaon energy calculated from the continuum dispersion relation 
and from directly htting the nonzero momentum two-point correlation functions on the 
ensemble discussed above. We do not observe any statistically significant deviations from 
the continuum dispersion relation. Further, while the statistical errors grow with increasing 
momentum, the kaon energies are consistent with a continuum dispersion relation within a 
2% statistical accuracy even at our largest simulated lattice kaon momentum. Therefore, 
we use the continuum dispersion relation to obtain the kaon energies at nonzero lattice 
momenta because this yields smaller statistical errors than the direct £t. 

The meson propagators from consecutive gauge-held configurations are, in principle, cor¬ 
related, so we look for possible autocorrelations by studying the effect of the block size on 
our fit results. We perform this test on every ensemble. As illustrated in Figure H] for two 
of the ensembles, the central values and errors are stable with increasing the block size, so 
we do not block the data or inhate the statistical errors in our analysis. 


12 


















1.852 
1.848 
1.844 
^ 1.84 

1.836 
1.832 
1.828 

123456789 

^min 


mass 
p-value 


1.0 

0.8 ^ 

^ 0.6 > 

a, 


- 0.4 

- 0.2 
n 




0.366 
0.3658 
0.3656 
0.3654 
0.3652 
0.365 

3 4 5 6 7 8 9 10 11 12 

^min 



FIG. 2. il-meson (upper) and kaon (lower) mass vs. tmin on the a ~ 0.12 fm , = O.lm)^ ensemble 

for fixed fmax = 15 and 30, respectively. The left and right vertical axes show the fitted mass and 
the p value (confidence level) of the fit, respectively. The filled circles show the values of tmin 
selected for the analysis. 


B. Extracting form factors 

We extract the lattice form factors /y, /x, and from the ratio of three-point to two- 
point correlation fnnctions. The three-point fnnctions are defined as 

T; fc) = e"“ »(Ox(0,0) y) OpT, x)), (3.4) 

T-,k)^Y. ''(Ojr(0,0) y) OpT, x)), (3.5) 

x,y 

where the kaon sonrce is at time slice 0 and the 5-meson sink is at time slice T. The 
sonrce-sink separations T are given in Table IVIIl Becanse we calcnlate the form factors in 
the 5-meson rest frame, only the kaon has nonzero momentnm k. 
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FIG. 3. + k'^) vs. kaon momentum in units of 27r/L on the a ~ 0.12 fm , = O.lm^ 

ensemble. The continuum dispersion relation is well respected through momentum 27r(l, 1,1)/T. 
The dashed lines show a power-counting estimate for the size of the momentum-dependent dis¬ 
cretization error for comparison. 


By inserting two complete sets of states, the three-point correlation fnnction Gg can be 
decomposed into snms over energy levels as 

k) = (3.6) 

m,n 


where 


-^mn 


{Q\Ok\K^^'^) 

2EP 




(E(^)|C>b|0 ) 

2mP 


(3.7) 


The contribntions from the first few terms dominate Gg at times snfficiently far from both 
the sonrce and sink. A similar decomposition applies to Gg^. 

We nse the averages introdnced in Ref. [l[ to snppress the contribntion from oscillating 
states in correlation fnnctions. We average the valne of the two-point correlator on snccessive 
time slices: 


G2(f+ 2) ' 

g-M™p-r2)_ 

(3.8) 


where Zp = (0|Gp|P) is the gronnd-state amplitnde of the kaon or B meson, and AMp is the 
energy difference between the gronnd and hrst oscillating state. For three-point fnnctions. 




( 0 ), 


C2{t) ^ 


C2{t) 2C2{t+l) 


+ 


Z'r} 


+ 


"p -M); 


2M 


( 0 ) 


+ G(AM^), 
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FIG. 4. Mb (left column) and Mk (right column) vs. block size on the a ~ 0.12 fm (top row) 
and a ~ 0.06 fm (bottom row), m[ = O.lm)^ ensemble. The fit results are stable as the block size 
increases. 


we also average the value of the correlator for two neighboring sink locations T and T + 1: 






+ + 1,T + l;fc) 


X 


Ct\t,T;k) ^ Cf''Ht,T+l;k) 


{T-t) g-E™ (t)g-M™ (T+l-t) 


g-E™(t+l)g-M®(T-t-l) 


+ 2, T; fc) + 2, T + 1; fc) 

H-^-- + 




g_E(^) (t+2) g-M4°’ (T-t-2) (t+2) {T-t-l) 


0{AEl, AEkAMb, AM] 


B) 


(3,9) 


(3.10) 
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We then form the ratios 




{t, T; k) ^ I ZE® 
k)C^{T-t) V e-^Kt ’ 


(3.11) 


where and are obtained from hts to Eq. fl3.2p with = \J+ k?. From 

Eqs. (I3.8p and (I3.10p . the ratio contains a f-independent term proportional to the 

desired matrix element, and other higher-order terms from the excited states. 

We show an example of the ratio on the a ~ 0.12 fm , = 0.1m( ensemble in 

Fig. 13 There is a short plateau region in the middle between 0 and T, with kaon excited- 
state contributions visible on the left and i?-meson excited-state contributions visible on 
the right. The i?-meson excited-state contributions, however, are smaller as indicated by 
the less dramatic falloff of the correlator on the right-hand side. We therefore choose to £t 
the correlator closer to the S-meson side including the contribution from a single S-meson 
excited state, but sufficiently far from the kaon that we can neglect kaon excited states. The 
£t function is given by: 

(t, T; k) = [l - (3.12) 




Q , Di, and AMb are £t parameters. 


where D 

models all excited states, we expect AMb 
excited state. 


Although the second term in Eq. fl3.12p 
to be close to the mass difference of the hrst 


1031 . Il04l | to Eq. fl3.12p . with priors determined 


We employ a correlated, constrained £t 
as follows. For the prior on we select a point from the middle of the plateau region 

and use its central value with the error inflated by a factor of two. For Hi, we use a prior 
of central value zero and width one. For AM^, we use the central value and width of 


- M 


( 0 ) 


— ivi^ obt aine d from the corresponding two-point correlator £t. 


augmented |104 


Xaug X 


E 


(ph) _ p{i)y 


a: 


We minimize the 


(3.13) 


where is the ith fit parameter, and P^*^ and cxj are the prior central value and width. 
We measure the goodness of fit using the y^^g/dof or p value, obtaining p from xlug 
the number of degrees of freedom equal to the sum of the number of data points and prior 
constraints minus the number of fit parameters. We choose the fit interval [tmin,^max] such 
that we obtain a good p value, using the same £t range for all momenta on the same ensemble. 
We select approximately the same physical £t ranges on the ensembles with different lattice 
spacings. Figure [5] shows sample hts of the three form-factor ratios on the a ~ 0.12 fm , 
m[ = 0.1m( ensemble. Figure [6] shows an example of the stability of the fit result against 
the variations of the fit range. We choose the preferred £t range to be [tmin,^max] = [8,12], 
where we find a good p value. The £t ranges and source-sink separations used on other 
ensembles are given in Table IVIIl 

To study the effects of residual excited-state contamination, we generated three-point 
correlators on the a ~ 0.12 fm, m[ = 0.14m'^ ensemble with several source-sink separations 
T = 18,19, 20, 21. We repeat the correlator hts with three sink-location combinations (T, T-|- 
1) = (18,19), (19,20), and (20,21), and the results are shown in Fig. [7] for four diherent 
momenta. We hnd no statistically signihcant diherences for all operators and momenta 
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TABLE VII. Pairs of source-sink separations T, T -|- 1 and fit ranges used in the fits. 


PS a (fm) 

r,r-M 

[^mim ^max] 

0.12 

18, 19 

[8, 12] 

0.09 

25, 26 

[10, 16] 

0.06 

36, 37 

[16, 24] 

0.045 

48, 49 

[20, 32] 


except for f± and /r at p = 27r(l, 0, 0)/L. These differences, however, are still snfficiently 
small that increasing the error on all p = 27r(l,0,0)/L points in the chiral-continunm £t 
does not change the physical form-factor resnlts. 

The £t parameters are proportional to the matrix elements {K^\J\B^). The lattice 
form factors are obtained as 


flat/p \ 

J± [Ek) - 

Pat/^ \ ^ MB + MKDo{k) 

^ ^/2M^ 


(3.14) 

(3.15) 

(3.16) 


The factor {Mb + Mk)/\/2Mb in fx in Eq. fl3.16p . which stems from Eq. fl2.7p . is evaluated 
with the physical meson masses to avoid introducing dependence not captured in the 
yPT formula. 


C. 6-quark mass correction 


The 6-quark hopping parameter used in our simulations differs slightly from the physi¬ 
cal value Kh because our production runs started before a more precise tuning of the 6-quark 
hopping parameter ki, was completed. For our desired accuracy, we need to apply a correc¬ 
tion. To this end, we have carried out runs with multiple values of on the a ~ 0.12 fm 
ensemble with = 0.2. In addition to the production value of = 0.0901, we repeated 

the run with = 0.0820 and 0.0860, allowing us to bracket the physical value Kb = 0.0868. 
The form factors depend on the 6-quark kinetic mass m 2 - At the tree level 


1 _ 2 1 
m2a mQa(2 -|- tuqu) 1 -|- rn'^a ’ 


(3.17) 


where 



(3.18) 


The values of Uq and Kcrit are given in Table IVIIII Following Ref. , we expand the form 
factor in about m' 2 . 


f{m'2,EK) = f{m2,EK) 


clln / 
clln m 2 




(3.19) 
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FIG. 5. Form-factor ratio fits on the a ~ 0.12 fm , m; = O.lm^ ensemble. From top 

to bottom, the three plots show the ratios for the temporal vector, spatial vector, and tensor 
currents. In the top plot, the data sets correspond to lattice kaon momenta k = 27r(0, 0,0)/L, 
27r(l, 0,0)/L, 27r(l,l,0)/L and 27r(l, 1,1)/L; nonzero momentum is required to extract the form 
factors in the bottom two plots, so there are only three sets of data in each of them. The gray 
horizontal bands show the fit results with statistical errors for in Eq. (13.121) . The black solid 

and dashed curves show the fit result within and extended beyond the fit range, respectively. 
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FIG. 6. Fit results of from different fit ranges on the a « 0.12 fm , = O.lm'^ ensemble 

with lattice kaon momentum p = ^(1,1,1). From top to bottom, the three plots show the ratios 
for the temporal vector, spatial vector, and tensor currents. We vary the fit range by changing tmin 
and tmax- The blue data point denotes the result from the fit range used in this paper. 
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FIG. 7. Form-factor ratio fits on the a ^ 0.12 fm, = 0.007/0.05 ensemble. From top 

to bottom, the three plots show the ratios for the temporal vector /y, spatial vector f±, and tensor 
Jt currents. The fit results for different pairs of source-sink separations T, T -|- 1 are shown as a 
function of Ek- The results from larger sink combinations are slightly displaced to the right for 
clarity. 
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TABLE VIII. The simulation and physical Hb We also include Kcrit and uq from the plaquette 
in this table for convenience, because they are used in the calculation of the 6-quark kinetic mass. 


Rs a (fm) 

am'Jam'^ 

i^'b 

Kb 

^crit 

Uq 

0.12 

0.01/0.05 

0.0901 

0.0868(9)(3) 

0.14091 

0.8677 

0.12 

0.007/0.05 

0.0901 

0.0868(9)(3) 

0.14095 

0.8678 

0.12 

0.005/0.05 

0.0901 

0.0868(9)(3) 

0.14096 

0.8678 

0.09 

0.0062/0.031 

0.0979 

0.0967(7)(3) 

0.139119 

0.8782 

0.09 

0.00465/0.031 

0.0977 

0.0966(7)(3) 

0.139134 

0.8781 

0.09 

0.0031/0.031 

0.0976 

0.0965(7)(3) 

0.139173 

0.8779 

0.09 

0.00155/0.031 

0.0976 

0.0964(7)(3) 

0.139190 

0.877805 

0.06 

0.0036/0.018 

0.1052 

0.1052(5)(2) 

0.137632 

0.88788 

0.06 

0.0018/0.018 

0.1052 

0.1050(5)(2) 

0.137678 

0.88764 

0.045 

0.0028/0.014 

0.1143 

0.1116(3)(2) 

0.136640 

0.89511 


where m 2 denotes the physical 6-quark kinetic mass. We determine the slope, in our 

companion work on the semileptonic decay B ttIu [^. Because the slope depends mildly 
on the daughter-quark mass, and the daughter-quark mass is tuned close to its physical 
value in our calculation, we neglect the daughter-quark dependence of the slope in this 
work. Finally, we quote of h, f±, and fr at the simulated daughter-quark mass as 

0.115(9), 0.139(13), and 0.126(13) j^. We hnd relative shifts due to b-quark mass tuning of 
about 0.5%-1.5% on the different ensembles. 


D. Chiral-continuum extrapolations 


The lattice form factors are computed numerically on ensembles with degenerate up- 
/down-quark masses that are heavier than the value in nature, as well as at nonzero lattice 
spacing. To obtain physical results, we hrst compute the form factors on several lattice 
spacings with varying up/down-quark masses and close-to-physical strange-quark masses, 
and then extrapolate to the physical light-quark mass and continuum (and interpolate to 
the physical strange -qua r k m ass) using heavy meson rooted staggered chiral perturbation 
theory (HMrSyPT) |l^.[Ti^. 

For the chiral-continuum extrapolation we use an HMrSyPT formula valid to leading 
order in l/rub and next-to-leading order (NLO) in the light-quark masses, kaon energy, and 
lattice spacing, supplemented by n ext-t o-next-to leading order (NNLO) analytical terms. We 
have tested both SU(3) HMrSyPT |l06l |. which includes the effects of dynamical pions, kaons, 
and T] mesons, and SU(2) HMrSyPT, in which the mesons with strange quarks are integrated 
out. In addition, we also consider hard-kaon HMrSyPT, which applies to semileptonic 
decays with energetic kaons. We hnd that NLO SU(3) HMrSyPT even supplemented with 
NNLO analytical terms, does not provide a good description of the data for /y lil-laj, 
and the p value of the ht is 10“®. On the other hand, SU(2) HMrSyPT describes the data 
well even at NLO. We therefore choose SU(2) HMrSyPT to perform the chiral-continuum 
extrapolations. 

The kaon energies in our numerical simulations are much larger than the rest mass of the 
physical kaon. Therefore standard HMrSyPT, which is derived for the situation in which 
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the kaon momenta are soft, may not provide a good description of our data throughout 
the available kinematic range. We therefore also consider hard-kaon HMrSyPT, which ap¬ 
plies for semileptonic decays with energetic kaons. Recently, Bijnens and Jemos derived 
the co ntinuum N LO hard-kaon (pion) HMyPT formulae for both B ^ K and B ^ tt pro¬ 
cesses [io^.lTn^. We derive the corresponding NLO staggered SU(2) and SU(3) hard kaon 
(pion) HMrSyPT formulae in Appendix It turns out that the chiral logarithms in NLO 
hard-kaon SU(2) HMrSxPT are identical to those in standard soft-kaon SU(2) HMrSxPT 
for B ^ K decays. This is likely the reason that the stan dard NLO SU(2) expressions 
describe our data even at such large kaon energies. Reference 109| found that the hard-pion 
theory can break down at three-loop level, but we only work at one-loop level here. 

The NLO SU(2) HMrSxPT formulae for R —)■ iP decays take the form 


d/2 


f\\ = 


(1 + logs) + + C\^\e + + C\^’xe 


hi). 


k2) 


(3), 


.(4), 




.- 1/2 


/± = 


Un{EK + 

(1 -|- logs) -|- C\_'^Xi + + C^I'^Xe + C[_'^Xa^ + 


UriiEx + ABj)ri 


(3.20) 


(3.21) 


where “logs” denotes nonanalytic functions of the light-quark mass and lattice spacing; the 
explicit expressions are given in Eqs. flA28D . flA33p . and flA34p . The dimensionless expansion 


P are 


2fimi 

“ SvrVf 

(3.22) 

2fimh 

M - g^2J2’ 

(3.23) 

a^A 

“ 87r2/r 

(3.24) 

\/2Ek 

XE = . r , 

4vr/^ 

(3.25) 


where A is the averaged taste-symmetry breaking parameter, a^A = A and /i 

denotes the leading-order QCD LEG; see Eqs. flA3D - flA14D for the dehnition. If HMrSyPT 
gives a good description of the data, we expect the / > 0, to be of order unity. The 
SU(2) yPT formulae do not contain explicitly; however, the low-energy constants (LECs) 
depend on rrih. Because the strange-quark masses on different ensembles are slightly different 
from each other, we include a term proportional to Xh in the set of analytic terms to account 
for the leading strange-quark mass dependence of the LECs and enable an interpolation to 
the physical strange-quark mass. 

Equations fl3.20p and fl3.2ip each contain a pole in Ek- The poles appear at negative 
energy — with 
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Ml 




(0) 




Mi 


2M, 


M. 


B’ 


B 


( 0 ) 


M, 


B- 


(3.26) 


The pole arises from low-lying states with flavor content hs and quantum numbers that 
depend upon the form factor: for and Jt, the relevant B* meson has = 1“, while for 
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TABLE IX. Fixed parameters used 

in the chiral fit [f 

p is the leading-order low-energy constant 

in QCD. and 

T\cP‘byij^ are the taste splittings and hairpin parameters for asqtad staggered 

fermions. 







a 0.12 fm 

a PS 0.09 fm a 

« 0.06 fm 

a « 0.045 fm 

continuum 

n/i 

6.831904 

6.638563 

6.486649 

6.417427 

6.015349 

r^a^Ap(10“^) 

0 

0 

0 

0 

0 

Aa(10 “^) 

22.70460 

7.469220 

2.634800 

1.040930 

0 

rfa^ Ap(10“^) 

36.61620 

12.37760 

4.297780 

1.697920 

0 

rfa^Ay(10“^) 

48.02591 

15.93220 

5.743780 

2.269190 

0 

A5(10“^) 

60.08212 

22.06520 

7.038790 

2.780810 

0 

r\a^5v 

0.00 

0.00 

0.00 

0.00 

0 


-0.28 

-0.09 

-0.03 

-0.01 

0 


/||, the B*q state has = 0^. In the chiral-continuum hts, we hx Mb to its experimentally 
measured value 5.27958 GeV (recall that we tuned the lattice 6-quark mass using the 
experimental i?s-meson mass.). We also use the experimentally measured value of the lowest- 
lying vector meson Mb* = 5.4154 GeV j^, which is stable apart from B* —)• for the 

pole position in the hts of f±_ and fr to Eq. fl3.21l) . Although a scalar B*q state has not 
been observed in experiments, theoretic al predictions estimate its mass to be just below 
the B-K production threshold [llO . 111 |. Therefore, in the ht of /y, we u se th e prediction 
Mb* = 5.711(23) GeV from a recent three-havor lattice-QGD calculation 53 for the pole 
position in Eq. fl3.2Up . 

Following the approach of Refs. 


1031. Il041 |. we constrain the parameters of the chiral- 


continuum ht with Bayesian priors and minimize the augmented xing dehned in Eq fl3.13l) . 
The chiral logarithms in Eqs. fl3.20p and fl3.2ip depend upon the universal B-B*-tt coupling 
which we constrain with a Gaussian prior of ce ntral valu e 0.45 and width 0.08. This 


prior is consistent with a direct l attice ca lculation |ll3 - 115 |. yet conservative enough to 


accommodate other lattice results liol. 117 |. The chiral logarithms also depend on the mass 


splittings between mesons of diherent tastes and on the leading-order LEG /r. These param¬ 
eters depend only on the light-quark action, and we hx them to the values determined in 
the MILG light-pseudoscalar analysis j3 ; see Table IIXI In the f\\ chiral-continuum extrap¬ 
olation, we account for the uncertainty on the scalar B* mass by taking a generous prior 
width of three times the theoretical error reported in Ref. 


or ±69 MeV. 


We constrain the coefficients of the LO and NLO analytic terms using priors 

with central values zero and widths two. To allow for higher-order contributions in the chiral 
expansion, we also include the complete set of NNLO analytic terms. These are proportional 
to xh XiXaX XiXe, XiXe: Xa?XE, Xa^Xs^ Xe: Xe: and We use prior central values of 0 
with widths 1 for the coefficients of the NNLO analytic terms. The systematic error from 
truncating the chiral expansion will be discussed in Sec lIVl 

Staggered yPT incorporates taste-breaking discretization effects from the light valence 
and sea quarks, but the lattice data also contain generic light-quark and gluon discretization 
effects as well as discretization effects from the heavy quark. We account for generic light- 
quark and gluon discretization errors by adding the term 2;as(aAQCD)^ in the HMrSyPT 
formulae with coefficient prior central value zero and width one. Similarly, to account for 
heavy-quark discretization effects in both the action and heavy-light currents, we add terms 
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of order and a^a with coefficients constrained by heavy-qnark power connting j^. At this 
order there are hve fnnctions (/s, /y, /s, /b, fx) that depend npon the bare heavy-qnark 
mass; their explicit forms are given in Appendix A of Ref. [I[. Dimensional analysis can be 
used to estimate the heavy-qnark error 

errorj oc /i(moa)(aA)'^™‘^'“^, (3.27) 


where /* is related to the mismatch between coefficients of the continuum operators in the 
action and currents and their lattice counterparts, and A is a typical QCD scale for heavy- 
light mesons which we take to be A = 500 MeV. As in Ref. , we add terms Zi x errorj to 
the HMrSxPT formulae for /||,±,t- The priors on the Zi have central values zero and widths 
equal to the square root of the number of times each function appears. (See Appendix A 
of Ref. [l|.) Because the discretization errors are included via the constrained fit in the 
chiral-continuum extrapolations, our results for the extrapolated form factors include the 
systematic uncertainties from light and heavy discretization effects. 

In summary, we use expressions derived in SU(2) HMrSyPT for the central chiral- 
continuum extrapolations of the form factors f\\, f±, and /t; these are shown in Fig. |8l 
The SU(2) theory describes our data well: the p values of the hts are 0.91, 0.94, and 0.98 
for /||, /j_, and /y, respectively. Our £t results for g.„ are 0.47(5), 0.46(4), and 0.47(3), re¬ 
spectively. At this stage, only the statistical, g.,r, chiral truncation, and discretization errors 
have been included. In the next section, we estimate the size of the remaining uncertainties 
before employing the 2 ; expansion in Sec. |V] to extend our results over the full kinematic 
range. 


IV. FORM-FACTOR ERROR BUDGET 

In this section, we estimate the systematic errors in the form factors, discussing each 
source of uncertainty in a separate subsection. We hrst discuss the error from the chiral- 
continuum extrapolation, which also includes heavy-qnark, light-quark, and gluon discretiza¬ 
tion errors. We then discuss the remaining systematic uncertainties from the heavy-light 
current renormalization, lattice-scale determination, light- and strange-quark mass determi¬ 
nations, hnite-volume effects, and 6-quark mass determination, discussing each in a separate 
subsection. As discussed previously, the systematic errors from g^, and heavy- and light- 
quark discretization effects are included in the statistical errors of the chiral-continuum 
extrapolation result through the constrained £t. Finally, we visually summarize the error 
budgets for the three form factors as a function of in Fig. [TT] 


A. Chiral-continuum extrapolation 

We use NLO SU(2) HMrSxPT supplemented by all possible NNLO analytic terms, as 
well as heavy-qnark, light-quark, and gluon discretization terms, in our preferred chiral 
extrapolations of f±, f\\, and /t. 

First, to estimate truncation effects, we compare £t results using NLO HMrSyPT, our 
preferred £t function with NNLO analytic terms, and the same £t function with the addition 
of the complete set of NNNLO analytic terms in Fig. |9l We see that the errors in the preferred 
£t with NNLO analytic terms are already saturated, since they are the same as the errors in 
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tiEk 



tiEk 


a ~ 0.12 fm — 0.20 

a ~ 0.12 fm — 0.15 

a ~ 0.12 fm — 0.10 

a ~ 0.09 fm = 0.20 

a ~ 0.09 fm — 0.15 

a ~ 0.09 fm = 0.10 

a ~ 0.09 fm — 0.05 

a ~ 0.06 fm — 0.20 

a ~ 0.06 fm = 0.10 

a ~ 0.045fm = 0.20 

continuum 



FIG. 8. Chiral-continuum extrapolations of /|| (upper left), fj_ (lower left), and /t (lower right) 
using NLO SU(2) HMrSxPT plus NNLO analytical terms. The squares, circles, triangles, and 
diamonds denote the = 0.2, 0.14, 0.1, and 0.05 data, respectively. The colored fit lines 

correspond to the different lattice spacings as indicated in the legend. The cyan band shows the 
continuum extrapolated curve with statistical error, which includes the systematic uncertainties 
due to Qn, and the heavy-quark, light-quark, and gluon discretization errors. Fit lines should pass 
through the data points of the corresponding color. 


the fit with NNNLO analytic terms. Hence, truncation effects are included in the statistical 
ht errors from our preferred ht. 


In addition, we also consider two alternative ht Ansatze for the chiral-continuum ex¬ 
trapolation. First, we consider NLO SU(3) hard-pion HMrSyPT, which provides a good 
description of our data, although the standard NLO SU(3) expressions do not. We use the 
result from the SU(2) HMrSyPT ht as our prefe rred ht, because the SU(2) theory converges 
faster than the SU(3) theory as studied in Ref. [n^ . We compare the ht results from NLO 
hard-kaon SU(3) HMrSyPT plus NNLO analytical terms and our preferred ht, and hud 
diherences between the central values of about 1-2% for all form factors and Second, we 


25 





















17 18 


19 


20 21 22 23 17 18 


(GeV^) 



(GeV^) 


19 20 21 

(GeV^) 


22 23 


FIG. 9. Chiral-continuum extrapolations with NLO, NNLO, or NNNLO analytic terms for /+ 
(upper left), /o (upper right), and fx (lower panel). In each plot, the grey band shows the statistical 
error from the preferred NNLO SU(2) The red and blue lines show the error from the fits 

with NLO and NNNLO analytic terms, respectively. 


consider the effect of the range of the lattice-QCD data to the extrapolated continuum 
result by omitting the k = 27r(l, 1,1)/L data from our fit. We find the differences are below 
1-2%. Figure [TUI summarizes the differences between the form factors obtained from the al¬ 
ternative chiral-continuum hts and the central results. Overall, the shifts of the continuum 
form-factor central values are within the quoted statistical errors of the preferred chiral £t 
that includes truncation effects. 


B. Heavy-light current renormalization 

To obtain the continuum form factors, we multiply the lattice form factors by the renor¬ 
malization constant given in Eq. fl2.17p . using the values of pj, Zy4^, and Zy4 listed in 

Tables MMl The statistical error on Zy^ is about 0.2%. By using the jackknife blocks 

^ss 
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FIG. 10. Deviations of alternate chiral-continuum extrapolations from the central results for /+ 
(upper left), /o (lower left), and fr (lower right). In each plot, the black curve shows the statistical 
error from the preferred NNLO SU(2) xPT. The blue and pink lines show the % difference from 
the central fit obtained by using SU(3) hard-kaon yPT and omitting k = 27r(l,l,l)/L data, 
respectively. 


of calculated on the same ensembles, we incorporate the statistical error from Zl /4 

bb Vfjb 

automatically in our £t results. 

The pj are calculated at one-loop order in perturbation theory. They are close to unity 
by design, since they are dehned as ratios of renormalization factors. Indeed their one-loop 
corrections are small, as shown in Table El We estimate the error due to truncating the 
perturbative expansion as in order to avoid sensitivity due to accidental cancella¬ 
tions. We obtain follows. For the scale-independent vector currents (1^* and V^), 

we simply look for the largest value of the one-loop coefficients for both currents on all of 
the ensembles. We find that the spatial vector current has a larger one-loop coefficient with 
Pvmax = O-l- evaluate Og at the a ~ 0.06 fm lattice spacing (the next to hnest), which 
yields an error of 1% for both components of the vector current. For the scale-dependent 
tensor current the perturbative corrections include logarithmic contributions due to their 
anomalous dimension, which are responsible for the growth of pT towards smaller lattice 


27 

























spacings seen in Table |V] . In order to estimate the truncation error, we remove the effect of 
the anomalous dimension by setting p = 2/a. We find that Pr^ax ~ which corresponds 
to a truncation error of 2% on px- In summary, we assign a perturbative truncation error 
of 1% on /+, 0 and an error of 2% on fx- 


C. Scale uncertainty 


We use Ti = 0.3117(22) fm in the continuum from Ref. [89 


to convert lattice quantities to 
physical units, where the quoted error includes both statistics and systematics. We repeat 
our analysis varying ri by plus and minus one standard deviation from its central value and 
use the larger change of each form factor as an estimate of the systematic error due to the 
scale uncertainty. We hnd differences of less than 1% for /y, f±, and fr throughout the 
region. 


simulated 


D. Light- and strange-quark mass uncertainties 


After the chiral-continuum £t, we evaluate the form factors at the physical quark masses 
rim = 0.000965(33) and = 0.0265(8) determined from the analysis of the light pseu¬ 


doscalar meson spectrum [4, l52|. We vary the quark masses by plus and minus one standard 


deviation and hnd the differences in all three form factors due to changing mi and m^ to be 
below 0.6% in the simulated region. 


E. Finite-volume effects 


The lattices used in this work have finite spatial volumes with > 4. We estimate 
the size of hnite-volume effects using HMrSyPT. In chiral perturbation theory, hnite vol¬ 
ume con tributions change loop-momentum integrals to sums which have been calculated 
in Refs. |l06 . 119 |. We employ continuum integrals in the preferred chiral-continuum ex¬ 


trapolations. To estimate the size of hnite-volume ehects, we evaluate the form factors 
with the LECs we obtain from the preferred chiral hts, and compare the results from the 
inhnite-volume formulae and the hnite-volume formulae on all ensembles used in this work. 
We try both SU(2) HMrSyPT and SU(3) hard-kaon HMrSyPT. We hnd that in all cases 
hnite-volume ehects are below 0.001%. Therefore, we neglect hnite-volume ehects in the 
total error budget. 


F. 6-quark mass correction 

We correct the form factors from the simulated to the physical Kb before we perform the 
chiral-continuum extrapolation. Including these corrections accounts for the dominant ehect 
from 6-quark mistuning, but small errors in the form factors remain due to the uncertainties 
in the K^-correction factors. The statistical errors in the slopes are at most about 10% 

for /_L,T at 27r(l, 1, l)/L, while the sizes of the Kb shifts applied to the data points are about 
l%-2%. We therefore take the systematic error from the Kb correction to be 2% x 10% = 
0.2%, which is conservative enough to accommodate the largest possible error in the shift. 
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FIG. 11. Statistical and systematic error contributions to /+ (upper left), /o (lower left), and /t 
( lower right). The left vertical axis label shows the squares of the errors added in quadrature, while 
the right vertical axis label shows the errors themselves. The filled, stacked curves from bottom to 
top show the total error when we add each individual source of error in quadrature one by one. 


G. Summary of the systematic error budget 

Figure [TTl visually summarizes the results for the statistical and systematic errors. For all 
three form factors, the combined chiral-continuum extrapolation error is the largest source 
of systematic uncertainty. The total errors in the form factors /+, /o, and Jt are below 5% 
for all > 17 GeV^, and are ~ 3% near g^ax- We quote numerical results for the form 
factors including all systematic errors over the entire g^ range in the following section, after 
the g^-extrapolation to the full kinematic range using the expansion. 


V. z EXPANSION OF FORM FACTORS 

The form factors obtained from the chiral-continuum £t are reliable for high momentum 
transfer, g^ > 17 GeV^. We only simulated kaons with momenta up to 27r(l, 1, l)/L, because, 
at higher momenta, the two- and three-point correlators become noisier and are subject to 
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larger discretization errors. Further, the HMrSxPT formalism used to take the continuum 
limit does not apply when Ek is too large. In particular, for Ek ^1-2 GeV the expansion 
parameter xe ^ G so the terms analytic in x'h increase with higher powers of n. Because 
of these limitations, a way to extend the form factors to high kaon energy, or, equivalently. 


= 0, is needed. In this paper, we follow Ref. |3l| and map to a new variable z such 
that 1^1 < 1. Constraints from unitarity, analyticity, and heavy-quark physics ensure that 
the expansion of the form factors in terms of z converges. Thus we can use the z expansion 
to obtain a model-independent parameterization of our form factors valid over the entir e 
kinematic range. This technique is now standard for analyzing B —)■ ttIu decays 


We hrst dehne the new variable 2 ; via the conformal mapping [31 






y/t+ — to 




(5.1) 


where t± = {Mb ± M^Y to is a free parameter that can be chosen to minimize \z 


for 
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the semileptonic-decay region. In this work, we use to = {Mb + Mk){\/Mb — \/MkY 
which maps the physical semileptonic decay region 0 < < 22.8 GeV^ to \z\ < 0.15. The 

small range of | 2 ;| helps control the truncation error in the 2 : exp ansion 
Using the new variable z, we expand the form factors as 
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U{q^) 

fo{q^) 

Mq^) 


P+{q 
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Po{q^: 
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m=0 
K-\ 

PT{q^) 

^ ^ ^ m=0 


T 


-(-1) 


m-K^ K 

K 


(5.2) 

(5.3) 

(5.4) 


The function P^ o,T{q^) = 1 — q^/M‘^ accounts for poles below and near the B-K production 
threshold. For the 2 ; hts of /+ and /t, we £x the location of the vector R* pole to the 
measured value Mb* = 5.4154 GeV |^. For the /o £t, we hx the loc ation of the scalar R*q 
pole to the lattice-QCD prediction Mb^q = 5.711 GeV from Ref. |ll2 |. We find that varying 
its location by three times the quoted theoretical error (±69 MeV) does not change the 
extrapolated form factor. 

The expression for /+ in Eq. fl5.2p was derived by Bourrely, Gaprini and Lellouch in 
Ref. [ 3 ^, and is commonly called the BGL parameterization. In the BGL expression for /_|_ 
in Eq. fl5.2p . the coefhcient of the term proportional to z^ is related to that of the lower- 
order terms. This constraint is due to the conservation of momentum and the analyticity of 
the form factors jssj. There is no analogous constraint for /q. We use the same expression 
for Jt as for /+ because they are proportional to each other at leading order in the heavy- 
quark expansion. These expressions were also used to analyze the lattice form factors for 
B —>■ KPl~ in Refs. [HI, 22 


Unitarity constrains the coefficients of the expansion such that 


V 

m,n=0 




(5.5) 
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where the values of Bmn are calculated using the Taylor expansion of the function discussed 
in Ref. 3^ and given in Table 1X1 We employ the same coefficients Bmn for /t and /+. The 


outer functio n d* defined in Ref. jl22| is used in the derivation of the Bmn- Although the 0 


of /o in Ref. (l22 was derived without a scalar pole, its form is not a ltere d by the presence 
of the pole, because \z\ always equals 1 on the unit circle. In Ref. 122l |. Becher and Hill 
showed that, in the limit of large 6-quark mass, the sizes of the z coefficients for /+ are even 
smaller tha n the expectation from fIS.Sp . Heavy-quark effective theory provides an estimate 


of the sum 122 




mn^m^n 

TT 


dt 


't+ 


t-to 


Tm 


— t 




(5.6) 


where i = -|-, 0, or T, and the 0 is an outer function. To calculate the integral in Eq. fl5.6p . 
we need to know the form factors in the range [t+, cx)]. For /_|_, we assume that f± gives the 
dominant contribution and has only the single B* pole. Taking the limit Mb —)■ C )0 gives 
the following simple form for /+(g^): 




Mb 

\/2 Mb 


f±{EK) 


Mb 

\/2Mb fwiEK + 


(5.7) 


We then use our determination of C® from our preferred chiral-continuum fit to obtain the 
estimate 


E BmnbmK « 0.012 . 


(5.8) 


This result means that Eq. fl5.5p is o nly a loose bound for /_|_. In addition, it is consistent 
with a power-counting estimate 122| |. which anticipates ^ BmnbmK to be of order (A/m;,)^. 
The analogous calculation for Jt gives a similar result. The analysis below will show that 
the heavy-quark (HQ) constraint on /_|_ (and Jt), Eq. 05.71) . together with the kinematic 
constraint, /o(0) = /+(0), suffices to keep the ^ fit under control. 

We assume a log-normal distribution on ^ Bmnbmbn to ensure that ^ Bmnbmbn is always 


positive. The contribution from this prior to the augmented is: 


xl. 


[In (^ ) Bmnbmbn^ P-] 




mn L'm 


(J^ 


(5.9) 


where p is the central value and a is the width of the prior. For /_|_ and /r, we choose p and 
a in Eq. 05.91) as ln(0.02) and ln(^^). This choice is conservative enough to accommodate 
the uncertainties in the estimates. 

We first generate from the continuum, physical quark-mass limit of the chiral extrapo¬ 
lation a few synthetic data points in the energy range of the simulated lattice data (g^ > 
16.8 GeV^). With the lattice spacing set to zero and the quark masses fixed to their physical 
values in Eqs. (I3.20p - (l3.2ip . the physical form factors depend upon at most six independent 
functions of the kaon energy E^- These are proportional to 1/{Ek + A^*), E^, Ek, E]^, 
E^, and Ej^. To the degree that the coefficients in front of these functions are correlated, 
the number of independent modes may be even fewer than six. If we generate too many 
synthetic data points, the covariance matrix will be singular. We therefore generate four 
synthetic data points each for /+, /o, and /r at = (22.86, 21.13,19.17,17.09) GeV^. These 
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TABLE X. Lowest-order coefficients B^n for B —)• Kl^l~ decay using Mb = 5.27958 GeV, = 
0.497614 GeV, and to = (-^s + Mk){VMb — The outer function used in the calculation 

is from Ref. [3^ with Xf+ = 5.025 x 10“^ and Xfo — 1-4:575 x 10“^. Although these y^s are derived 
for the B —)• ttIu process, the calculation in Ref. shows the difference between y^s of the 
B —)• Kl^l~ and B —)• nli' process is less than 10%. Therefore, we quote the inputs from Ref. 321 
to obtain these Bmn- All Bmn not listed here can be obtained from the relations Rm(m+n) = Bon 
and Bran — Bnm- 



Boo 

Bqi 

Bo2 

0 

CO 

Ho4 

Bob 

f+,T 

0.0161 

-0.0003 

-0.0104 

0.0002 

0.0022 

0.0002 

h 

0.0921 

0.0132 

-0.0483 

-0.0168 

- 0.0001 

0.0024 


cover the simulated lattice-momentum range and are approximately evenly spaced in q‘^. We 
also fit with synthetic data from a smaller and larger range and hnd consistent results. 

The full covariance matrix of the synthetic data points includes both the statistical and 
systematic error: 

/ofull /ostat I /osyst fK ^ C\\ 

^mn ^mn ^mn ’ 1^0.iUJ 

where m, n denote the four values. The systematic error contribution is calculated as 


^syst \ ^ i i 

^mn / j ^ n 


(5.11) 


where the index i runs over the sources of systematic error discussed in Sec. HVl Because 
we assume that the systematic errors are 100 % correlated between values, all nontrivial 
correlations between points are due to statistical fluctuations of the chiral-continuum £t 
results. 

We first fit /+, /q, and fx simultaneously in a combined fit using K = 3 (three free 
parameters) in Eqs. fl5.2p - fl5.4p without any constraints on the coefficients. Table IXll presents 
the results of these fits. We plot the fit results in Fig. [121 Although we do not impose the 
kinematic condition /+(g^ = 0 ) = /o(g^ = 0 ), it is approximately satisfied with separate 
hts. Adding HQ constraints on the /+ and Jt fit makes the results even more consistent 
with the kinematic condition (see Fig. [12]), and reduces the errors on /_|_, at low g^. We 
then fit /+, /o, and Jt simultaneously with the kinematic constraint, and still including the 
HQ constraints on /_|_ and Jt, which further decreases the extrapolation error in the form 
factors at low g^. We implement the kinematic constraint by setting a prior of /+ — /o at 
g ^=0 with central value zero and width of 0 . 00001 . 

We show the Bmnbmbn bootstrap distribution of /+ and fr from two fits with and 
without the HQ constraint in Fig. [131 Adding the HQ constraint moves the distribution 
of ^ Bmnbmbn to Smaller values. We also compare the ^ Bmnbmbn distribution of /o from 
two fits in Fig. [TTl One is a fit with /o only, the other is a combined /+ and /o fit with 
the kinematic constraint. Adding the kinematic constraint decreases ^ B, 
separate /o fit. Again, the result shows that the unitary constraint on the ^ B, 
is a loose bound. 

We also check the truncation error by repeating the fit with K = 4. Because in K = 3 
fits, the coefficients 62 are not well determined by data, and the results are zero within error, 


mnbmbn from the 


-'mnbmbn 


of /o 


we add a prior of 0(2) on b\ coefficients as in Ref. [2]J to control the fluctuations of the 


higher-order terms. All of the coefficients b from fits with K = 3 and 4 are summarized 
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2 ; 2 

FIG. 12. Separate z-expansion fits of /+, /o (left) and /r (right) without (upper) and with 
(lower) HQ constraints on the sum of coefficients for /+ and Jt- The synthetic data points are 
generated at large (small z) in the region of simulated lattice momenta. The kinematic condition 
f+{q‘^ = 0) = /o(g^ = 0) is satisfied better when the HQ constraint is applied to /+. (Recall that 
the factor = 1 at = 0.) 


in Table IXIi The results from different K are consistent with each other. The coefficients 
are zero within error and have little impact on the central value of the hnal result. We 
therefore conclude that the ^ truncation error is well controlled. 

We record our hnal, preferred results from K = ?> z hts including both the heavy-quark 
and kinematic constraints in the third column of Table IXII and we give the corresponding 
correlation matrix in Table IXIIl Together with the pole masses (also in Table IXip and 
Eqs. fl5.2p - fl5.4l) . this information allows the reader to reconstruct our form-factor results 
throughout the full kinematic range. Our hnal form-factor results as a function of 2 ; and 
are plotted in Figs. fTSHT^ 
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FIG. 13. Histogram of the the sum of coefficients Bmnbmbn for /+ and /t from fits with and without 
the HQ constraint. Use of the HQ constraint moves the distribution of Bmnbmbn to smaller values. 



FIG. 14. Histogram of the sum of coefficients ^ Bmnbmbn for /o from an independent fit and from 
a combined fit with /+ that imposes the kinematic constraint at = 0. 
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TABLE XL Results of z-expansion fits of the B ^ K form factors /+ (top panel), /o (middle panel), 
and fr (lower panel) using the formulae defined in Eqs. (|5.2I1 - (I5.4I) with to = {Mb + Mk){VMb — 
^/M^f 0, Mbi = 5.4154 GeV in /+,t, Mb*^ = 5.711 GeV in /o, Mb= 5.27958 GeV, and Mk = 
0.497614 GeV ff. 



unconstrained 

constrained 



HQ 

HQ -|- kinematic 


K = 3 

K = 3 

K = 3 

A = 4 

K 

0.437(22) 

0.451(20) 

0.466(14) 

0.466(15) 

K 

-1.41(33) 

-1.15(27) 

-0.89(13) 

-0.89(16) 

bt 

-2.5(1.4) 

-1.4(1.1) 

-0.21(55) 

-0.19(61) 

bt 

- 

- 

- 

0.3(1.1) 

Bmnbmbn 

0.16 

0.07 

0.02 

0.03 

/+(0) 

0.18(10) 

0.256(80) 

0.335(36) 

0.336(44) 


0.285(11) 

0.286(11) 

0.292(10) 

0.292(11) 

6? 

0.19(14) 

0.20(13) 

0.28(12) 

0.28(13) 

6° 

-0.17(49) 

-0.15(48) 

0.15(44) 

0.18(68) 


- 

- 

- 

0 .2(1.7) 

Bmnbmbn 

0.02 

0.02 

0.02 

0.02 

MO) 

0.309(39) 

0.311(38) 

0.335(36) 

0.336(44) 

bl 

0.440(25) 

0.453(23) 

0.460(19) 

0.459(20) 

bT 

-1.47(37) 

-1.17(30) 

-1.09(24) 

-1.11(24) 

bl 

-2.7(1.6) 

-1.4(1.2) 

-1.11(97) 

-1.15(95) 

bl 

- 

- 

- 

-0.2(1.1) 

y^ Bmnbmbn 

0.18 

0.07 

0.05 

0.05 

/r(0) 

0.17(11) 

0.254(87) 

0.279(67) 

0.276(68) 

p value 

0.57 

0.39 

0.34 

0.97 


VI. TESTS OF QCD PREDICTIONS FOR FORM-FACTOR RATIOS 


Because lattice-QCD calculations of the B ^ K semileptonic form factors have until 
recently been unavailable, theoretical calculations of S —)■ Kl^l~ observables sometimes use 
expectations from heavy-quark symmetries to relate them to others that can be constrained 


from experiment or computed with QCD models (see, e.g. Ref. 123| ). Heavy-quark symme¬ 


try is also commonly used in pheno menologi cal calculations of the related decays B —)■ ttZ+Z” 
B — )■ and B — )■ K*'y [H, 123l - 129l |. Here we use our lattice-QCD form factors to 

directly test these heavy-quark symmetry relations in R —)■ A decay at both high and low 




35 






















FIG. 15. /+, /o, and Jt z-expansion fits. The synthetic data points are generated at large (small 
z) from LECs of the HMrSyPT fit result. The kinematic constraint f+{q'^ = 0) = /o(g^ = 0) is 
applied exactly in the combined /+ and /o z-expansion fit. The vertical dashed lines correspond 
to q‘^=0. We use three coefficients [IT = 3 in Eqs. (I5.2p - (|5.4I) ] for /+, /o, and fr- 



EIG. 16. /+, /o, and /r vs. q^ based on the z expansion. The kinematic constraint /+(g^ = 0) = 
foiq'^ = 0) is applied exactly in the fit. We use three coefficients [K = 3 in Eqs. (I5.2p ~ (l5.2p ] for /+, 
/o, and /t. 


A. Low-recoil predictions from heavy-quark symmetry 


In the soft-kaon {Ek -C Mb) and chiral limits, the vector and scalar fo rm factor s can be 
related using heavy-quark effective theory and chiral perturbation theory |l3n . 131 1: 


/o 

hm = 
q2^Ml /+ 


IBs 

fsi 


q^M'- 


Bt 




( 6 . 1 ) 


where the decay-constant ratio accounts for heavy-quark corrections of 0{l/mb)- Heavy- 


quark spin symmetry relates the vector and tensor form factors in the soft-kaon limit as |123 
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TABLE XII. The coefficients bi from the z-expansion fit (the first line) and their correlation matrix. 
The upper index +, 0, and T denote the form factors They are from the z-expansion fit 

formulae defined in Eqs. (I5.2I) - (I5.4I) . We use to = {Mb + Mk){-\/Mb — [s^], M^ = 

5.4154 GeV in /+,T; Mb*^ = 5.711 GeV in /o, Mb= 5.27958 GeV and Mk = 0.497614 GeV p. 



b^ 

bt 

b} 

OO 

6? 

bl 

b^ 

bf 

bf 

Mean 

0.466 

-0.885 

-0.213 

0.292 

0.281 

0.150 

0.460 

-1.089 

-1.114 

error 

0.014 

0.128 

0.548 

0.010 

0.125 

0.441 

0.019 

0.236 

0.971 

K 

1 

0.450 

0.190 

0.857 

0.598 

0.531 

0.752 

0.229 

0.117 

K 


1 

0.677 

0.708 

0.958 

0.927 

0.227 

0.443 

0.287 

bt 



1 

0.595 

0.770 

0.819 

-0.023 

0.070 

0.196 

b^o 




1 

0.830 

0.766 

0.582 

0.237 

0.192 

6? 





1 

0.973 

0.324 

0.372 

0.272 

60 






1 

0.268 

0.332 

0.269 

bl 







1 

0.590 

0.515 

bf 








1 

0.897 

bl 









1 


Ill: 


= ^^{^A 

q^^Ml U 


Mb{Mb + Mk) 

q2 


+ 0{A/mb ), 


( 6 . 2 ) 


where the scale-dep enden t coefficient ^(/i) incorporates cor recti ons of 0 {a‘^) to the leading 


Isgur-Wise relation |l33l | and is given in Eq. (2.5) of Ref. 123 


We can estimate the size 
of higher-order corrections in the heavy-quark expansion from power counting. Taking 
A = 500 MeV and rrib = 4.2 GeV gives A/rnt ~ 12% and (A/m^)^ ~ 1%. Equations fl6.ip 
and fl6.2l) also receive corrections from the kaon recoil energy that are of 0{Ek/TU b). For 
Qmax ^ ^ 14 GeV^, this ratio varies from 12% < Ex/mb < 40%, so such corrections are 

expected to be signihcant even at low kaon recoil. 

Figure [T71 left, compares the quantity (%//+) x (1 —obtained from our lattice 
form factors with the theoretical prediction Eq. fl6.ip . For the theoretical estimate , we take 


fs*/fss = 0.953(23) from the recent four-flavor lattice-QGD determination in Ref. jl34l | and 
g.^ = 0.45(8) as in our chiral-continuum £t. The width of the theoretical band is from the 
uncertainty on and does not include any other errors. Figure [T71 right, compares the 
quantity {fx/f+) x {q^)/{M b{Mb + Mk)) obtained from our lattice form factors with the 
theoretical prediction Eq. 06.21) using rrib = 4.18 GeV and as^{mb) = 0.2268, such that 
K{mb) ~ 0.88 |l23 . 124j |. We do not show any errors on the theoretical prediction. 

The observed lattice form-factor ratios /o//+ and /t//+ at are lower than the 
theoretical expectations by 38% and 15%, respectively; by = 14.5 GeV^ the differences 
grow to 51% and 46%, respectively. Although the observed disagreement with the theoretical 
expectation for the tensor form-factor ratio is large, it is within the size expected (from 
simple power counting) for higher-order corrections due to the kaon recoil energy. The 
scalar form-factor ratio, however, differs from the theoretical expectation by a much larger 
amount. In Fig. 25 of Ref. we compare the quantity (%//+) x (1 — for the 

related decay R —)■ ttZz/ with the heavy-quark prediction in the soft-pion limit. The observed 
agreement near q^^^ is better, which suggests that the discrepancy is indeed due to the light 
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FIG. 17. Comparison of lattice form-factor ratios with theoretical predictions from heavy-quark 
symmetry at low recoil. Left: (/o//+)/(l — q^/Mp , vct sus from lattice QCD (red curve 
with error band) and heavy-quark symmetry plus yPT jl3Cll | (gray horizontal band). The width of 
the theoretical band includes the uncertainty on = 0.45(8) but no other theory errors. Right: 
{It/ f+) X {q‘^)/{Mb{Mb + Mk)) versus from lattice QCD (red curve with error band) and the 
improved Isgur-Wise relation 123l ] (black horizontal line). 


pseudoscalar-meson recoil energy, which is larger for B ^ K than for i? —)• vr. Thus our 
lattice form-factor results suggest that one should be cautious in using heavy-quark relations 
derived in the soft-pion/kaon limit for phenomenological predictions, especially for decays 
with K or K* hnal-state mesons. 


B. Large-recoil predictions from QCD factorization 


In the large-recoil limit {E^ S> Tf^), heavy-quark sy mme try relates the vector, scalar. 


and tensor form factors to a single universal form factor jl35 


1 . fo 

hm — 

Ek'>Mk J-i- 

lim ^ 

Ek^Mk /-)- 


2E 


K 


Mb 
Mb ■ 


-b 0{A/mb), 

- Mk 


M, 


+ (9(A/mQ. 


(6.3) 

(6.4) 


B 


The 0{as) corrections to th ese l eading large-recoil expressions were derived using QCD 
factorization (QCDF) in Ref. 136| . and the resulting expressions are given in Eqs. (62)-(63) 
of that work. Higher-order corrections in the heavy-quark expansion are expect ed to be 
about A/mb ~ 12%, while corrections to the QCDF predictions from Ref. |I^ are 

expected to be about 5%. 

Figure [18] comp ares the lattice-form-factor ratios with the theoretical large-recoil pre¬ 
dictions from Ref. jl36| . For the O(aA corrections, we take the decay constants = 
190.5(4.2) MeV from FLAG [l2l| and fx = 156.2(7) MeV from the PDG We 
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FIG. 18. Comparison of lattice form-factor ratios with theoretical predictions from heavy-quark 
symmetry at large recoil. Left: (/o//+) x Mb/{2Ek) versus from lattice QCD (red curve with 
error band) and theoretical prediction with 0{as) corrections |l36l | (gray curve with error band). 
Right: {/t/ f+) X {Mb)/{Mb + Mk) versus fro m lattice QCD (red curve with error band) and 
theoretical prediction with 0{as) corrections 136| (gray curve with error band). 


take the first inverse moment of t he 5-meson distribution amplitude A^^(2.2 GeV) = 
[0.51(12) GeV]-^ from LGSR [mJ, where the quoted theory error covers the spread of 
other determi nations from QGD/light-cone sum rules and the operator-product expan¬ 
sion |l38l - 14n |. We take the hrst and second moments of the kaon distribution ampli¬ 
tude ai{2 GeV) = 0. 061(4) and af^(2 GeV) = 0.18(7) from a recent three-flavor lattice- 
QGD calculation |l4lj . We use our own determination of f+{q‘^ = 0) = 0.335(36). We 
take as^{mb) = 0.2268 as described above and q:s^( 2.2 GeV) = 0.279 |l42| . The left 
panel of Fig. [TSl shows the quantity (/o//+) x {Mb)/{2Ek), while the right panel shows 
(/t//+) X {Mb)/{Mb + M^)- The widths of the theoretical bands in Fig. [18] are from the 
uncertainty on and f+{q^ = 0), and do not include any other errors. 


For {fo/fj^)xMB/{2 Ek), the lattice-QGD result differs from the theoretical predictions by 
at most 1%, which is well within the expected size of heavy-quark corrections. For {Jt/ f+) x 
{Mb)/{Mb + M^), the lattice-QGD result is marginally consistent with the theoretical 


expectation of Ref. jl36 


A more recent NNLO calculation within soft-collinear effective 

The new 


143 


theory updates the large-recoil predictions to include 0 {a‘^) corrections 

= 0 result for fx/ f+MB/{MB+MK) = 0.817 is in better agreement with the ratio obtained 
from lattice QGD. Overall, the uncertainty on the lattice-QGD tensor form-factor ratio at 
low is too large to draw any quantitative conclusions. (The vector and tensor form factors 
are not strongly correlated at low g^.) Thus, while the scalar form-factor ratio suggests that 
the large-recoil predictions may be reliable, some caution is nevertheless warranted in their 
use for phenomenology given the limited number of tests they have undergone. 
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VII. SUMMARY AND OUTLOOK 


As discussed in Sec. El Table IXlTI presents our final results for the form factors /+(g^), 
/o(g^), and for the semileptonic process B —)■ Kl^l~. These entries, which consist 

of the coefficients of the BCL expansion, Eqs. fl5.2p - fl5.4p . together with the correlations 
among them, can be used to reconstruct our form factors with errors for all values of 0 < 
< ?max- This information can also be used to compute form-factor ratios and (differential) 
rates with squares of linear combinations of the form factors. 

Figure [11] shows a comparison of our results with others in the literature. A t = 0, 
our result is consistent with a light-cone-sum-rule result from Khodjamirian et al. |l44| . For 
all our results are consistent with the only other unquenched lattice-QCD calculation 
from the HPQCD Collaboration j2l|. Our form factors are somewhat more precise than 
HPQCD’s, especially at high g^, because we used more ensembles with finer lattice spacings 
and lighter quark masses. The total errors, including both statistical and systematic errors, 
are less than 4% at high g^, and at low g^ about 10% for /+ and 30% for fx- 

More generally, our results can be used to compute any B —)■ Kll observable, including 
asymmetries and decay rates, for all possible dilepton final states (Z = £, r, z/), and even 


lepton-flavor-violating modes |145| . We present a tho roug h analysis of observables for B —)■ 


K semileptonic decays in a companion publication |l46j |. where we also present ratios of 
observables for B —)■ Kll to i? —)■ nil decay processes. The three form factors /+, /o, and fx 
suffice to parameterize the factorizable hadronic contributions to B ^ K semileptonic decays 
in any extension of the Standard Model. Other hadronic uncertainties, such as violations 
of quark-hadron duality due to intermediate charmonium resonances, must, of course, also 
be reliably estimated to obtain complete Standard-Model and new-physics predictions for 
B ^ K processes. If deviations from the Standard Model are observed in any B —)■ Kll 
decay channel, accurate results for the form factors will be essential to disentangling the 
underlying physics. 

The main sources of uncertainty in our form factors are from the chiral-continuum ex¬ 
trapolation and extrapolation to low g^. We plan to reduce these uncertainties with newer 




FIG. 19. Our form factors compared with light-cone-sum-rule results 144l | and the other un¬ 
quenched lattice-QCD calculation [21]. 
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gauge-field ensembles that are being generated by the MILC Collaboration |l47l 148 1. These 


ensembles use the highly-improved staggered quark (HISQ) action for the light, strange, and 
charm quarks. This action is designed to have small er di scretization effects which will help 


reduce the size of the continuum extrapolation errors (l49 |. In addition, the HISQ ensembles 


include ensembles with physical pion masses, which will eliminate the need for the chiral 
extrapolation and the associated errors. Indeed, these ensem bles have already been used 

and charmed-meson |l5l| physics. In particular, 
that an analysis with physical and 


to improve the precision for kaon |l50 


we found with D- and H^-meson decay constants |151 


unphysical quark masses provides better statistical precision with no penalty in systematic 
errors. 
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Appendix A: S —Kll form factors in SU(2) S^PT 


We use expressions derived in HMrSxPT |l06 | as the low-energy effective theory of QCD 
in which the degre es of f reedom a re pions and kaons for our chiral-continuum extrapolations. 


SU(3) HMrSyPT jl06|, Il52l. Il53| was applied to i? —)■ Trlu semileptonic decays Q . More 


recently, SU(2) HMyPT |ll8l. Il54ll55| was also considered as an alternative effective theory 
in studies of heavy meson physics. We derive the SU(2) HMrSyPT formulae for form factors 
calculated with staggered quarks in this appendix. These formulae can be used for B ^ irlu, 
B —)■ Kl^l~ and D-meson semileptonic decays. Our r esult s are consistent with earlier studies 
of HMyPT for continuum QCD and Wilson quarks |l55 | after taking the continuum limit 
of the HMrSxPT expressions. The differences in the detailed expressions can be absorbed 
into redehnitions of the scale or LECs. 


1. /|| and /j_ in SU(3) HMrSxPT 


The SU(3) HMrSyPT expression of /y is the same as in Eq. fl3.2Up . We only list the 
expression of the log terms for unitary points /j = ^ = mi^s) here. 

For the B ^ n process, the chiral logs in SU(3) HMrSxPT are given by |l06 |: 




logS||,SU(3) 




1 - 39 = 


[2Ji(m^=) + hirnK,E)] 


(47r/)2 [ 16 V L 2 

-h 2/2(771^ =, E) + E) 

1 -|- “^9^ j f N 1 r / N 

H-1- 

'KqI-i) 


-hirrijy) - 2 l 2 {mjy, E) 


+ [E ^ A] 


(Al) 
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For the B ^ K process, the chiral logs in SU(3) HMrSxPT are given by: 


R^K 


loS®||,SU(3) 


(4„/)2 I 16 - 3 dli(.;m,.s) + i/i(ms,E) 


E 


+ 2/2(mx,s, E) + l2{ms,E, E) 

- ]^Ii{ms,i) + + hirunj, E) - himsj, E) 


a^'y 


him^'y) - h{m^y) + him^'y, E) - him^y, E) 


"^ri'y "^vy 


- ^ Rf’^^ {{msy, m^y, {m^y}) f ^Ii{mjy) + hirnjy, E)j 

4 cSQrn-r,n ^ / 

2 Rf^H{^ny,mrfy,mr,>y}]{msy}) hirrijy) 


je{S,r),ri'} 

^91 

ie{7r,r?,r7'} 

+ [1/ ^ A] 


off-diagonal meson masses are: 






(A2) 

appear 

in Eqs. 

flAlD and flA2D are as follows. 

The flavor 

= 

fi{mi + 

mi) 

+ a^As, 

(A3) 

mly = 

fi{mi 

mg) 

-|- ci^A=, 

(A4) 

rnl,E = 

fi{mg + 

mg] 

1 -|- o^A=, 

(AS) 


where mi and are sea-quark masses and the taste label E has values P, V, T, A and I. 
The masses of flavor-neutral mesons in the taste vector channel are 


mloy = mfjy = mjjy = n{mi + mi) + a^Ay, 
m^gy = jji{ms + mg) -|- a^Ay, 

^Iv = \ 

^ f 2 2 3 

I m.ir ,, -I- m A ,, -I- 

2 


"^li'y — 9 ('^uy + ^%y + + E j , 




^iv - 


x2 2 

y) - —{m-^^,-m 


2 \ I 9(o^^\, 
uyj 


vJ 


(A6) 

(A7) 

(AS) 

(A9) 

(AlO) 


2 V-5,y ■-u,v/ ' 

The taste-axial case just requires substituting A for V. For the taste-singlet case, we have: 

= ^u,i = ^D,i = + ^i) + (All) 

mil = ^{mh + mh) + a^A/, (A12) 

A 

'' 3,1 


2 , 2^2 

T^n,i = + 

ml,j = ml. 


(A13) 

(A14) 
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The momentum integrals Ji and I 2 that appear in the chiral log terms are defined as: 


hijn) 
him, A) 


Fix) 




+2A2, 


\/l — tanh ^ (\/l — xh 0 < x < 1, 
— \/x^ — ltan“^ (\/x2 — 1) X > 1, 


(A15) 

(A16) 

(A17) 


where A is the renormalization scale. _ 

Similarly, fj_ on the unitary points in NLO SU(3) HMrSyPT is given by 


f± = 


+ 


C(o) 

~T 


E + A% + D 


C(o) 


fiE + As* 


-(logs + (A18) 


where Ab* = Mb* — Mb- (The SU(3) expression has one extra chiral log term D comparing 
with the SU(2) expression we used in our analysis.) There are two chiral log related terms 
par amet erized by D and logs in Eq. flAlSIl . For the i? —)■ tt process, the SU(3) expressions 
are 


^SV{3) 


'islE 


-t [2 Jr^m,,=, E) + Jr\mK,s, B)] 


(dvr/)^ 1 16 
-ljr\m^j,E) + ^jr\mr,,i,E) 

+ ^ i-a^6v)Rf’^^ i{m„y, m^y, {msy}) Jl'^^irujy, E) 

3&{'E,rix'} 


+ [V^ A] 


(A19) 


and 




logS_L^gU(3) 


ly 

(47r/)2 I 16 


1 + 3^^ 
2 

1 


[2Iiim,,y) + himKy)] 


E) + igpr'l™,,/. E) + 


Miim-^j) - /i(m^,/) 


y o?5yRf^\{m^y,mr,y,mr,'y]-,{msy]) 

jG{ 7 r,r), 77 '} 




+ IV^A] 


(A20) 
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For the B ^ K process, the SU(3) expressions are: 


nB^K 

^SV{3) 


+ i-a‘^5y)Rf^^ {{msy, m^y, m^^y}] {m^y}) E) 

ie{ 5 , 7 ?,T;'} 

+ [V^A\\, 


(A21) 


and 




^ \-y 

{4:71 fy I 16 


9 -I- 1 

- 2/i(ms,s) - 3glh{my 


E) + ^h(my - + l/,(ms,/) 

^7T I TSUb/^ 7Ti\ TSUb/ 


+ a"h 


2 x/ 
V 


^9l 


^I'y - ^Iv 

?[3,1] 


jrym,y,E)-jrym,>y,E) 


^ {{m^y, m^y, m^^y}; {msy}) him^y) 


j^R,v,v'} 


n ^f’^H{^sy,mrjy,m^>y}]{m^y}) Ii{mjy) 


j€{S,ri,r)'} 

+ |1/^ A]l . 


(A22) 


The dehnition of the meson mass terms and R are the same as for the /y case. The f± 
expression has an extra fnnction Ji, that is dehned as: 

Ji(m, A) = (^-m^ + In + ^(A^ - - yA^ + ^rri^ , (A23) 


J™ (m, A) = Ji(m, A) 


27rm^ 

3A 


(A24) 


2. /|| and fj_ in SU(2) HMrSxPT 

We derive the SU(2) formnla for /y and f± based on the SU(3) expression. We also nse 
the same expression for /t as for f± as discnssed in Sec. IIIIDl To obtain the SU(2) limit of 
an SU(3) expression, we treat the strange-qnark mass as inhnitely heavy. The SU(2) form 
does not contain explicitly, bnt all LECs depend implicitly on m^. Becanse onr lattice 
data have slightly different on different ensembles, we keep the analytic term which is 
proportional to m^. Next, we consider all terms in the SU(3) chiral log expression. If a term 
is proportional to or Inm^ in the large limit, it is absorbed into the redehnition of 
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other LECs. If a term is proportional to l/m* or 1/ Inm^ in the large nis limit, it does not 
appear in the SU(2) expression. We now derive the form of the chiral log terms in the SU(2) 

limit. 

For the chiral log terms in /y, because we take rUg to inhnity, all related terms, such 
as mK,E, iTirij and rrij^iyiAii, go to inhnity. They are absorbed into LECs. Only rUrfy/A 

is hnite and goes like \J We now consider all contributing chiral log terms. 

• Ii{m) goes like m^lnm^, so only /i(Mn-) survives. 

• him, E) diverges as 27imE when m —)■ oo, so only hiM^,, E) survives. 

• The ratio 

hjmr^iy) - Iiimr,y) + hjmr^'y, E) - hjmr^y, E) 

diverges as 2 In rus at large m^, so it is removed. 

• We hnd that 


lim 

ms—>-oo 


„2 r/ p[3,l] 
a OyiA-ttj 


i{msy, mrjy, mr,>y}; {m^y}) 


4, J = S 

-^ = 0, J=V ■ (A25) 

-4, j = r(, 


When this term multiplies Ji or I 2 , it is divergent as —)■ cx) for j = S or j = rj'. 

For j = Tj, the Ji and h are hnite, but the total contribution is zero as uig —)■ 00 . So 
these terms are removed. 


• We hnd that 


2 , 


lim i{m.^y,mr^y,mr^>y}-,{msy}) = “2, 


5 ;? 


J =71 

j = V , (A26) 


V/A _ 


Am'g 


= 0 , j = 7]' 


SO only j = 71 and j = rj terms contribute in the SU(2) theory. 

In summary, for the B ^ 71 process, the chiral log in SU(2) HMrSyPT is given by 


logs 


&’(2, = I S E [(1 - SglWrn^s) + 2/.(m, =, E)] + ,) 


3 ( 9 ; - 1) 


him ^ y ) - 2/2(771^, y,E) 


^igl - 1) 


himr,y) - 2/2(m^,y,E) 


+ [E ^ A] 


(A27) 
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For the B ^ K process, the chiral log in SU(2) HMrSxPT is given by 


logsK) = 

_l—^ [2/i(m,r,v) ~ 2/i(m^y)] + [1/ —)■ A] 


(A28) 


We then derive the expression for the f± chiral log terms in SU(2) HMrSyPT. We use 
the same treatment of analytic terms as was done for /y. To calculate the SU(2) chiral log 
terms, we consider the large irtg limit of Ji. 


lim Ji{m, E) —>■ —m^lnm^ —?• —oo, 


lim 

mg—>-oo 


Ji{mny,E) - Jl'^'^{mn',v,E) 


^l',v - ^Iv 


—)■ 2 Inm, —)■ cxD. 


(A29) 

(A30) 


So all Ji related terms are absorbed into the redehnition of LECs and disappear. 

Via a procedure similar to that for /y, we obtain the SU(2) chiral log terms in f± for the 
S —)■ TT channel: 


Di 


- [2jf‘(,n^y, E) - -B)] + |r ^ yl] I , 


(A31) 


^B—¥7r 


log®±,SU(2) 


E 


(47r/)2 I 16 ^ 
1 + 3^7^ 


1 + 377^ 


[2Ji(m^,s)] 




12 


3Ji(m^p) 


+ 


2 I E) + 


2 I E) + 




Similarly, the SU(2) chiral log terms in B ^ K are: 

dB^k _ r. 

-^SU(2) — U, 


(A32) 


(A33) 


logs!:^u(2) = 

_l—^ [2Ji(m7ry) — 2/i + [V —)■ A] 


(A34) 


Equations flA3ip . flA32p . and flA34p are written with a structure similar to their SU(3) 
counterparts, which makes it easier to implement a unihed computer code for the various 
choices of yPT studied in this paper. 
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3. Form factors in hard-pion/kaon ChPT 


The hard-kaon (pi on) conti nuum HM^PT for B ^ K and i? —)■ tt semileptonic decays 


was derived in Refs. [Tn7l.H 


The pion or kaon with large E is integrated out from the 
theory and its effects are absorbed into the LECs. We derive the hard-kaon (pion) limit of 
the HMrSxPT in this section. We hrst study the asymptotic behavior of the integrals which 
contain or E^- We hnd that 


2 

TTL 

him.E) ^ AoE‘^\n{E'^) + AiE'^ + A 2 \nE - m^\n{—) , (A35) 

E) BoE'^ ln(E2) + B^E"^ + R 2 InE + Rg , (A36) 


in the large E limit, where the coefficients Ai and Bi are either constants or analytic functions 
of m. The divergent terms in the large E limit decouple from the expression. The analytic 
terms in m are absorbed into the redefinition of the LECs. So the rules to derive the 
hard-kaon (pion) HMrSyPT are the following: 


Replace the term l 2 {m,E) by —/i(m) 
Remove E) term 


To compare our results with Refs. [I07l. Il08l |. we set all taste splitting parameters, hairpin 


parameters and lattice spacings to zero, 
(pion) HMyPT results. 


We then can reproduce the continuum hard-kaon 


iogs!:^u(3) = -(7 + - (i + '^gl) 


,3 9 2x/i(m^) 

A ' (47r/)2 

— a ' 

'4^ 


,1 3 

'2 ' 2 ' 


h{mK) 

(4vr/)2 


T2 


(4vr/)2 


1 9\-^ll^^7rj /-I 2 2\ 

logS±,SU(3) = -(^9n) (A^f\2 “ (1 + -jg-^) 


himK) 

{Airfy 


D = 0 . 




>-71 


[n + 75'i 




uum, which is also found in Refs. liH ITnsI 


1 ^ 

4^"^ (47r/)2 

, (A37) 

’ (4vr/)2 ’ 

(A38) 


(A39) 

= logsj^su(^) in 

the contin- 
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